Band insulator to Mott insulator transition in a bilayer Hubbard model 
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The ground state phase diagram of the half-filled repulsive Hubbard model in a bilayer is inves- 
tigated using cluster dynamical mean field theory. For weak to intermediate values of Coulomb 
repulsion U, the system undergoes a transition from a Mott insulating phase to a metallic phase 
and then onto a band insulating phase as the interlayer hopping is increased. In the strong coupling 
case, the model exhibits a direct crossover from a Mott insulating phase to a band insulating phase. 
These results are robust with respect to the presence or absence of magnetic order. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.10.Hf, 71.27. +a, 71.30.+h, 71.45.Lr 
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Strong Coulomb interaction can localize electrons even 
in the absence of disorder to give rise to a gapped elec- 
tronic state known as the Mott insulator. This state is 
of tremendous interest because upon addition of carriers 
it produces high temperature superconductivity and in- 
triguing features such as a pseudogap in the normal state 
of the cuprates. Band insulators, on the other hand, lack 
the ability to conduct because they have either a filled or 
an empty band. The topic of transitions between band 
and Mott insulators has received considerable attention 
recentlji^j^jiiJiiSi!, because it provides a useful tool to 
understand the fundamental differences between the two 
states. Recent progress in the control of Fermionic atoms 
trapped in optical lattices adds to the experimental in- 
terest as it would be feasible to investigate the evolution 
of correlated ground states in a many body system. In 
a recent study, it was shown that the ground state of a 
Mott-insulating Haldane phase in a generalized Hubbard 
ladder can be connected adiabatically to a band insulat- 
ing phase 8 . On the other hand, studies of the two di- 
mensional ionic Hubbard model (IHM), which contains 
an alternating local potential ±A on the A or B sub- 
lattice, find an intervening phase between the band and 
Mott insulating phases rather than a continuous evolu- 
tion when the Coulomb repulsion is varied. Cluster dy- 
namical mean field theory (CDMFT) calculations for the 
IHM at zero temperature suggest that the intervening 
phase is a bond-ordered insulator— whereas finite tem- 
perature Quantum Monte Carlo techniques propose it to 
be a metalS. 

Another theoretical model for the study of the band 
to Mott insulator crossover is the bilayer Hubbard model 
(BHM). The presence of copper-oxide bilayers, well de- 
scribed by doped Hubbard planes, in many of the high- 
T c materials adds to the practical importance of study- 
ing the bilayer Mott system theoretically. The BHM 
was studied recently using Dynamical mean field the- 
ory (DMFT),— and a smooth crossover between band 
and Mott insulators was observed. Since this study 
dealt with the infinite dimensional lattice, intralayer spa- 
tial correlations are absent. As a two-dimensional Hub- 
bard model exhibits fundamentally different behavior 
from the infinite-dimensional Hubbard model, such as 
a pseudogap 9 and d-wave superconductivity^ upon dop- 



ing, a study incorporating the two-dimensional nature 
of the Hubbard planes in the BHM is expected to show 
richer physics. 

In this work, we study the evolution of an antiferro- 
magnetic (AF) Mott insulator into a band insulator as 
the interlayer hopping is varied in a bilayer Hubbard 
model at half-filling using cluster dynamical mean field 
theory 11 . With interplane hopping set to tj_, the Hamil- 
tonian for the bilayer Hubbard model can be written as 

H = -t (cl aa Cj (TO , + H.c)j + U y n ila n^ a 

{ij)aa ia 
iera 

Here a = 0, 1 labels the two Hubbard planes and U rep- 
resents the onsite Coulomb repulsion. Ci aa denotes the 
destruction operator at site i, with spin a and on plane 
a. Chemical potential is chosen as fj, = U/2 so that the 
density is fixed to unity. In the limit U — ► 0, it can be 
seen that the system goes from an uncorrelated metal to 
a band insulator at t±_ = At when the splitting between 
the bonding and antibonding bands produced by the two 
layers results in a finite gap. At large U and t± — > 0, the 
system decouples into two Hubbard systems with a finite 
Mott gap and AF long range order in each layer. 

In the recent study on the infinite-dimensional BHM, 
the two-plane system was reduced to two impurities em- 
bedded in a self-consistently determined bath£. For 
t± = the authors found a phase transition at U = U c 
from a metallic phase to a Mott insulating phase. As 
t± is switched on, the Mott insulator (MI) evolves con- 
tinuously into a band insulator (BI) with no sign of a 
phase transition. In what follows, we shall argue that 
in-plane spatial correlations provide a competing energy 
scale that is missing in the infinite dimension. Taking 
both intraplane and interplane spatial correlations into 
account within CDMFT we find that for U < 8t as t± is 
increased from zero the MI goes through a metallic phase 
before entering a band insulating phase. For larger val- 
ues of U, the system undergoes a direct crossover from a 
MI to a BI without an intervening metallic phase. 

CDMFT is a non-perturbative technique where the 
full many-body problem is mapped onto local degrees 
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FIG. 1: Tiling the infinite bilayer with a 2 x 2 cluster 



of freedom treated exactly within a finite cluster that is 
embedded in a self-consistent bat h 11 ' 12 . It is a natural 
generalization of single-site DMFT— to incorporate spa- 
tial correlations. The method has passed rigorous tests 
for the ID Hubbard model where it compares well to 
exact solutions^ and has been applied to a variety of 
problems^. The first step in the CDMFT method in- 
volves a tiling of the infinite lattice by a finite cluster. 
As seen in Fig. [1] we consider a 2 x 2 tiling cluster which 
includes both the in-plane and intcrplanc hopping pro- 
cesses with a unit cell defined by the vectors u = x + y 
and v = x y. Using CDMFT, the BHM on the infinite 
lattice reduces to the cluster-bath Hamiltonian below, 
that is subject to a self-consistency condition: 



H.C.) + Uj2 n ^ 



E 



€ i n a & m (j & m a 
ma irtfia 



J2 V ^(o-l a c^ + H.c).(2) 



Here, fi, v = 1, • • •, N c denote indices labeling the cluster 
sites and m = 1, • • •, JV& represent those in the bath. The 
self-consistent calculation proceeds by an initial guess 
for the cluster-bath hybridization V mM(T and the bath 
site energies e mcr to obtain the cluster Green's function 
G{f . Applying the Dyson's equation, S c = Gq 1 — G~ l , 
where Go denotes the non-interacting Green's function 
for the Hamiltonian in Eq. @, the cluster self-energy is 
obtained. S c is then used in the self-consistency condi- 
tion below to determine the local Green's function for the 
lattice: 



Gloc(z) 
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dP u dP v 
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z + /i-f(P u ,P v ) -£<.(«) 

(3) 

Here, P u and P v denote reciprocal lattice vectors con- 
jugate to the spatial unit vectors u and v respectively, 
while z = iu) n is the Fermionic Mastubara frequency. 
Gioc can be used to obtain a new Weiss field, Go, and 
consequently, a new set of Vm^o and e mff 's by means of 
a conjugate gradient minimization program to close the 
iterative loop. In this work, the Lanczos method is used 
to solve the cluster-bath Hamiltonian and we fix the bath 
size to Nf, = 8. The Lanczos method can access both the 
strong and weak coupling regimes with equal ease and it 
is well suited to compute dynamical quantities directly 




FIG. 2: Phase diagram for bilayer Hubbard model. Shaded 
portion shows the extent of the metallic region (paramagnetic 
metal (PM) and antiferromagnetic metal (AFM)) in the pres- 
ence of magnetic order. Dashed line showing the magnetic 
phase boundary extends into the insulating region separating 
the paramagnetic Mott insulator (PMI) from the antiferro- 
magnetic Mott insulator (AFMI). Suppression of magnetic 
order extends the metallic phase down to lower t± shown by 
the dotted curve. 



in real frequency. Rotational symmetries of the cluster 
on a square lattice, together with particlc-holc symmetry 
at half-filling, reduce the number of bath parameters sig- 
nificantly. These bath parameters can depend on spin, 
allowing for symmetry breaking solutions. For details of 
the method, we refer to earlier wor k 11 ' 14 . 

Let us start from the phase diagram with suppression 
of the magnetic order. Numerical results are shown in 
Fig. We observe three states; a Mott insulating state 
at large U, a band insulating state at large t±, and a 
metallic state in between. At U = 0, the BI and the 
metallic states are separated at t± = 4 where the bond- 
ing and antibonding bands are separated. Critical t± 
decreases with increasing U because the width of bond- 
ing and antibonding bands broaden due to correlation. 
Above U = 12, the BI and the MI states are adiabati- 
cally connected, but by close inspection of the charge gap, 
two regions are distinguishable; in the MI (BI) phase, 
A c decreases (increases) with increasing t±. Those fea- 
tures are consistent with the DMFT result on the infinite- 
dimensional BHM. In contrast to the infinite-dimensional 
case, the MI state extends down to small U and t± in 
the two-dimensional case because of the intralayer spa- 
tial correlation, and the system passes through two phase 
transitions when increasing t± at < U < 12. Thus, the 
two-dimensional BHM provides topologically distinct be- 
havior as compared to the infinite-dimensional one. As 
will be discussed in detail below, this essential feature is 
not altered even by including the magnetic ordering. 

We next investigate the magnetic properties of the bi- 
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FIG. 3: Staggered magnetization a function of t± for U — 
At, 6t, 8t, I2t, 16t 



layer Hubbard model in detail. The staggered magne- 
tization M = 7if — Jij, is plotted as the interplane hop- 
ping is varied in Fig. [3] for different values of U. At 
t± = 0, the two planes are decoupled and the value of 
the staggered magnetization increases as U is increased, 
but quickly saturates to a value below unity because of 
quantum fluctuations. With increasing t±, M stays fairly 
constant until t± = t and then drops rapidly to zero at 
about t± = 2t. This value of t± = 2t is universal and 
does not depend on U for all U > 8t. The vanishing 
of the staggered magnetization at small U is caused by 
an increase in the effective bandwidth and a consequent 
reduction in the exchange coupling. In contrast, critical 
value t± w 2t at large U is due to the instability to- 
ward the formation of a local singlet. In this limit, the 
BHM is reduced to the 5 = 1/2 Heisenberg model with 
the intralayer exchange J = 4t 2 /U and the interlayer ex- 
change J± = At\/U. By taking a dimer on the interlayer 
bond as a unit and introducing a mean-field decoupling 
on the intralayer bonds, critical exchange coupling for 
the transition from a Neel ordered state to an interlayer 
dimerized state is obtained as J± = 4J at T — 0, thus, 
t± = 2t. Staggered magnetization M = 2(S Z ) at t± < 2t 
is also computed. As shown with a dashed line in Fig. [3j 
this mean-field approximation to the bilayer Heisenberg 
model gives reasonable agreement with CDMFT at large 
U. 

Next, we look at the local density of states for the BHM 
as a function of t± at U = 8t as shown in Fig. 2] As t± is 
increased the Mott gap reduces till it closes completely at 
t± = 2t to give rise to a metallic state. The LDOS shows 
a clear quasiparticle peak along with the lower and upper 
Hubbard bands characteristic of Mott physics at t± = 
2 At. Upon further increase of t± we see that a gap opens 
abruptly to reveal a band insulating state. Increasing 
t±_ in the band insulating phase causes the charge gap 
to grow monotonously unlike the Mott insulating phase 
where it reduces to zero. In Fig. [5] we plot the double 
occupancy and the charge gap as a function of t± for 
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FIG. 4: Local density of states as t±_ is varied at U = 8t. Spin- 
resolved density of states shown for t± = 1.0 where magnetic 
order exists. Note the metallic state at t± = 2.4t. 




FIG. 5: Double occupancy (top panel) and charge excitation 
gap (bottom) as functions of t± for various U indicated. 



different values of U. For U < 8t, as ij_ is increased 
in the Mott phase, the double occupancy rises and then 
shows a sharp kink as the system enters a metallic phase 
and also upon going into a band insulating phase. For 
larger values of U, as expected, the crossover from a Mott 
insulator to a band insulator is smooth with the double 
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occupancy showing no signs of a kink. The charge gap for 
large U shows a continuous decrease in the Mott phase 
but never goes to zero. The crossover to a band insulator 
around tj_ = 2t is easily identified by the increase in the 
charge gap with tj_- 

The phase diagram for the bilayer Hubbard model that 
emerges for our calculation is shown in Fig. [5] We iden- 
tify a region in the phase diagram which shows an inter- 
mediate metallic region, depicted by the shaded portion, 
which separates the band and Mott insulating phases. 
The width of the metallic region narrows with increasing 
U and eventually for U > 8t there is a direct crossover 
from a band insulator to a Mott insulator. The mag- 
netic phase boundaries are depicted by the dashed line 
which separates a antiferromagnetic metal from a para- 
magnetic metal. In the insulating portion of the phase 
diagram the magnetic phase boundary separates the an- 
tiferromagnetic Mott insulator (AFMI) from a param- 
agnetic Mott insulator (PMI). As we discussed earlier, 
emergence of intermediate metallic region separating BI 
and MI is a generic behavior in two-dimensional BHM 
independent of magnetic ordering. 

In summary, we have obtained the zero temperature 
phase diagram of the bilayer Hubbard model at half- 
filling using cluster dynamical mean field theory. We 
computed the double occupancy and the local density 
of states as a function of the interlayer coupling ij_ and 
the Coulomb repulsion U to identify a clear demarcation 
between a Mott insulating phase and a band insulating 



phase as the interlayer coupling is varied. For weak to 
intermediate values of interaction (U < 12t), the Mott 
insulator at small tj_ is separated from the band insula- 
tor at large t± by an intermediate metallic phase. For 
U > 12t, there is a direct crossover from a Mott insula- 
tor to a band insulator as t± is increased. This crossover 
is most clearly manifest in the behavior of the charge 
gap which first reduces in the Mott phase with increas- 
ing tj_ and then increases in the band insulating phase. 
These results should be contrasted with the DMFT re- 
sults for the infinite-dimensional model which always find 
a smooth crossover between the BI and MI states* 5 These 
differences arise due to the presence of in-plane spatial 
correlations in the two-dimensional model stabilizing a 
Mott phase for all U at half-filling which is accessible by 
cluster DMFT. We completed the phase diagram by fur- 
ther allowing for magnetic order. The phase diagram is 
essentially unchanged, except that the onset of the metal- 
lic phase occurs at a larger value of t± for a fixed value 
of U. 
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